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Abstract. We report on X-ray and optieal observations 
of a ROSAT X-ray source, RX J1605. 3+3249, selected 
from the all-sky survey on the basis of its spectral soft- 
ness and lack of bright optical counterpart. The ROSAT 
PSPC energy distribution is well fitted by a blackbody 
with kT = 92±6eV and Nh = l.l±0.4x lO^o cm-^. X-ray 
observations spanning 6.5 years fail to reveal any flux or 
spectral variability on any time scales. The ROSAT HRI 
error circle only contains a R = 23.25 M star which is un- 
likely to be associated with the X-ray source. We conclude 
that RX J1605. 3-1-3249 is a probable nearby isolated neu- 
tron star detected from its thermal emission. The present 
data do not allow to unambiguously determine the X-ray 
powering mechanism, cooling from a young neutron star 
or heating by accretion from the interstellar medium onto 
an old neutron star. However, the long term stability of 
the X-ray flux favours the young neutron star hypothesis. 

Key words: X-ray general. Stars: neutron. Stars: indi- 
vidual: RXJ1605.3-f3249 



1. Introduction 

Several arguments based on present metallicity of the in- 
terstellar medium, rate of supernovae and on properties 
of the observed radio pulsar population indicate that of 
the order of 10^ to 10^ old isolated neutron star (INS) 
should exist in the Galaxy. Ostriker, Rees & Silk (1970) 
were the first to propose that a sizeable fraction of these 
old and radio-quiet neutron stars could be heated by ac- 
cretion from interstellar medium and be again detectable 
by their far-UV and soft X-ray emission. Early modelling 
of this population by Treves & Colpi (1991) and Madau 



& Blaes (1994) led to the conclusion that re-heated old 



neutron stars should indeed appear in large numbers in 
soft and UV all-sky surveys with a possible concentration 
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in the directions of highest interstellar medium densities, 
namely the galactic plane in general and more specifically 
molecular clouds. 

Boosted by these predictions several optical identifi- 
cation campaigns of ROSAT X-ray and UV sources were 
initiated and it readily became clear that the number of 
possible isolated neutron star candidates was substantially 
below average model predictions (e.g., Manning et al. 



1996, Motch et al. 1997, Banner 1998). However, four good 



isolated neutron star candidates were discovered so far in 
the ROSAT all-sky survey. These candidates share as com- 
mon properties a soft X-ray spectrum with black body 
temperatures below 100 eV, no radio emission detected 
and very high Fx/Fopt ratio in excess of 10^. The X-ray 
brightest of these candidates are RXJ1856. 5-3754 which 
was optically identified with a V=25.6 blue object (Wal- 
ter et al. |1997| ) and the pulsating source RX J0720. 4-3125 
(Haberl et al. |1997 ) which has no counterpart brighter 
than B=26.1 (Motch & Haberl |l 994 Kulkarni & van Kerk- 
wijk |l998D Other good cases are RX J0806.4-4123 (Haberl 
et al7]l99^ ) and IRXS J130848.6-K212708 (Schwope et al. 
1999| ). 

In the mean time, the possibility that a fraction of 
these candidates could rather be young neutron stars has 
gained considerable credit. The lack of detectable radio 
emission from these sources could be explained by a B, 
Pspin position behind the radio pulsar death line or more 
simply by beaming effects which are stronger at long spin 
periods (Wang et al. 1998). Finally, RXTE observations 
have demonstrated that soft 7-ray repeaters are newly 
born neutron stars with extremely high magnetic fields 



(Kouveliotou et al. 1998) and belong to the population of 
magnetars proposed by Duncan & Thomson (1992). As ra- 
dio emission can be quenched by the strong magnetic field, 
these objects may remain undetected by classical radio 
means and their birth rate could amount to 10% of that 
of ordinary pulsars (Kouveliotou et al. 1994). Since the 



primary store of energy in a magnetar is that in the mag- 
netic field, B decay could constitute a significant source of 
heat, allowing magnetars to remain detectable in X-rays 
over longer times than ordinary pulsars (Heyl & Kulkarni 
19981) . 
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In fact, several of the INS found so far could be young 
neutron stars, perhaps descendant of soft 7-ray repeaters, 
rather than old accreting INS as originally thought. De- 
tailed study of the few known INS and determination of 
their X-ray powering mechanism is therefore of high im- 
portance. 

In this paper, we report on X-ray and optical observa- 
tions of one of the X-ray brightest isolated neutron star 
candidates. The source R X J160 5. 3-1-3249, also known as 
RBS 1556 (Schwope et al. 1999 ), was extracted from the 



ROSAT all-sky survey on the basis of its soft spectrum 
and lack of bright optical and radio counterpart. 

2. Selection of isolated neutron star candidates 
from the ROSAT all-sky survey 

Thanks to its soft sensitivity well suited to the detection 
of sources with Tbb = 20-100 eV, as expected from young 
cooling neutron stars or from old ones re-heated by accre- 
tion, the ROSAT all-sky survey offers a highly valuable 
database for detecting these elusive objects. In order to 
find candidates, we selected ROSAT all-sky survey sources 
displaying HRl and HR2 hardness ratios compatible with 
intrinsically soft spectra slightly modified by a reasonable 
amount of interstellar absorption. Hardness ratios 1 and 
2 are defined as 

(0.5 - 2.0) - (0.1 - 0.4) 



HRl = 



HR2 = 



(0.1 - 0.4) + (0.5 - 2.0) 

(1.0 -2.0) -(0.5 -1.0) 
(1.0-2.0) 



where (A-B) is the raw background corrected source count 
rate in the A— B energy range expressed in keV. Based 
on simulations of black body energy distributions folded 
with the ROSAT PSPC response we decided to extract 
all-sky survey sources having hardness ratios compatible 
(i.e. within one standard error value) with HRl < -0.25 
and HR2 < -0.5. This parameter space corresponds to Tbb 
< 100 eV and Nh < 4x10^1 to 3xl0^° cm-^ for Tbb = 
40 eV and Tbb — 100 eV respectively. Our observational 
strategy was then to identify all sources in large sky ar- 
eas down to the faintest X-ray flux level possible in order 
to find new candidates and also efficiently constrain the 
space density of these objects. Results from this global 
study will be presented in a later paper. The possible INS 
nature of RX J1605. 3-1-3249 = IRXS J160518.8-I-324907 
was discovered while identifying the northern sample. 
IRXS J160518. 8+324907 has a count rate of 0.875±0.041 
cts/s, HRl = -0.70±0.03 and HR2 = -0.58±0.10. To 
our knowledge, RXJ1605. 3+3249 is the brightest isolated 
neutron star candidate in the northern hemisphere. 



out during the 1998 PSPC revival period from 1998 Febru- 
ary 18 till 22 for a total exposure time of 4413 s. The sec- 
ond observation was performed with the HRl from 1998 
March 2 to 4 and lasted 19307 s. The source was detected 
on both occasions. 

3.1. X-ray spectral analysis 

The time averaged PSPC spectrum is well represented by 
a blackbody energy distribution with best fit parameters 
of Tbb = 92 eV and Nr = 1.1 10^° cm'^ {^li = 50.9). 
At the 95% confidence level, the allowed range is Tbb 



86 - 98 eV and Nr = 0.6 - 1.5 lO^'^cm- 



We show 



in Figs. |l| and ^ the best blackbody fit and the corre- 
sponding allowed spectral parameter range. The observed 
fiux corresponds to a bolometric luminosity of Lboi = 
l.lxlO'^^ (d/lOOpc)^ erg s~^. Assuming isotropic black- 
body emission, the source radius scales as i? = 1.1km 
(d/lOOpc). Black body temperature and line of sight ab- 
sorption compare well with those observed from other INS 
candidates such as RXJ1856. 5-3754 (Tbb = 57±le V, Nh 
= 1.4+0.1 10^0 cm-2; Waher, Wolk & Neuhauser |l996| ) 
and RXJ0720.4-3125 ( Tbb = 79±4eV, Nh = 1.3±0.3 



lO^o cm-2; Haberl et al. |1997| ) 



3.2. Search for X-ray variability 

Source flux appears remarkably constant over years to 
weeks time scales. We list in Table |l| the various source 
intensity measurements which are all compatible with 
a strictly constant emission. For the HRl observation 
we computed an equivalent PSPC count rate using the 
blackbody spectral description derived from the 1998 
pointed PSPC observation. The energy distribution mea- 
sured from hardness ratios did not change either between 
the almost 6.5 years elapsed from survey to pointed PSPC 
observations. 

We also searched for both aperiodic and periodic vari- 
abifity in the ROSAT PSPC and HRl time-series. In aU 
cases we failed to detect any signiflcant variability. Applied 
to light curves binned in 100 s intervals the Kolmogorov- 
Smirnov test gives a 95% confidence upper limit of 43% 
on variability amplitude for the PSPC time series and 
does not provide useful constraints for the HRl data. 
Power spectrum analysis also fails to detect any peri- 
odic signal with an estimated upper limit of 50% and 
36% full amplitude modulations for periods longer than 
1 s in the PSPC and HRl data respectively. We note that 
if RX J1605. 3+3249 had exhibited pulsations with an am- 
plitude simil ar to those seen in RX J0720.4-3125 (24% , 
Haberl et al. 1997), we would not have detected them. 



3. ROSAT observations 

ROSAT observed the field of RX J1605. 3+3249 in pointed 
mode on two occasions. The first observation was carried 



3.3. Source position 

Considering the extreme optical faintness of the counter- 
part any attempt to identify the source obviously requires 
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Table 1. Flux measurements 

Date Mode Detector PSPC cts/s Note 

10-11 Aug 1991 survey PSPC 0.875 ± 0.041 HRl = -0.70±0.03 HR2 = -0.58±0.10 

18-22 Feb 1998 pointed PSPC 0.900 ± 0.021 HRl = -0.66±0.02 HR2 = -0.68±0.04 

02-04 Mar 1998 pointed HRI 0.925 ± 0.020 HRl cut rate = 0.21 ± 0.01 cts/s 
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Fig. 1. Best blackbody fit to the PSPC count distribution 
of RXJ1605.3+3249 with Tbb = 92eV and Nh = 1.1 lO^o 
cm"^ 



as good as possible X-ray localization. Although the po- 
sitioning of the X-ray source on the HRI or PSPC instru- 
mental reference grid can be accurate at the arcsec level 
for relatively bright sources such as RX J1605. 3-1-3249, the 
uncertainty on the attitude of the satellite introduces a 
dominant 8-10 arcsec error. However, if enough identi- 
fied and well localized sources are present in the field of 
view, it is possible to correct for the unknown attitude 
error and retrieve the intrinsic accuracy achievable with 
the given detector. 

Above a Maximum Likelihood of 8, a total of 23 and 14 
sources are detected in the PSPC and HRI fields of view 
respectively. We have cross-correlated the HRI source list 
with SIMBAD, FIRST and USN0-A2 catalogues. Among 
these detections, 4 PSPC and 6 HRI sources have a posi- 
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Fig. 2. contour plot showing the allowed parameter 
space for Nh and Tbb at the 1 , 2 and 3 a confidence levels 



five match in the searched catalogues. We also took into 
account the 4" shift across the 40' HRI field of view due 
to the pixel siz e bein g 0.9972±0.0006 instead of 1 arcsec 
(Hasinger et al. |199^ ). For one of the HRI and PSPC field 
source, RX J1605. 5-1-3239, we had two possible identifica- 
tions, either the nucleus of the spiral galaxy CASG 1345 or 
the FIRST radio source located 6.5" away. Identifying the 
HRI source with the FIRST entry yields attitude correc- 
tion vectors incompatible with those derived from other 
sources in the HRI field of view. We therefore assumed 
that the HRI source was identified with the galactic nu- 
cleus. We fitted to the differences between X-ray and op- 
tical positions expressed in arcsec, relations of the form 



ax 



^opt 



5" (XiMA - Xccnter) (1 " 0.9972) + Cte 
and (5x -<5opt = 5" (Yjma - Ycentcr) (1 - 0.9972) + Cte 
where Xjma and Yima are the position of the sources on 
the grid of 5" size HRI pixels. The fit took into account 
the error on the positioning on the instrumental grid. We 
show in Fig. ^ the best fit obtained for the declination 
axis. 

Applying these relations to RX J1605.3-f 3249 moves 
the uncorrected HRI position by 0.4" to the East and 3.6" 
to the North to a = 16h 05m 18.66s and 6 = -1-32° 49' 
19.7" (2000.0 eq.) with a 1 cr error of 0.64". Using the 4 
identified PSPC sources yields a compatible position at a 
= 16h 05m 18.75s and S = +32° 49' 17.8" (2000.0 eq.) with 
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Table 2. BRI photometry of objects A,B,C and D derived 
from CFHT observations 



400 



Fig. 3. X-ray to optical position offset as a function of 
HRI pixel number 



a 1 CT error of 3.4". The HRI and PSPC pointed positions 
are well within the 90% confidence ROSAT survey error 
circle (a = 16h 05m 18.8s and S = +32° 49' 07.5" (2000.0 
eq.) with a. 1 a error of 7.0"). 

4. Optical data 

4-.1. Observations 

First optical observations took place from 1998 April 14 to 
17 using the Canada France Hawaii telescope. Images and 
spectra were obtained with the OSIS V instrument which 
provides image stabilisation and active guiding. With the 
2048 X 2048 STIS2 CCD the pixel size is 0.156" on the 
sky. We acquired several 15min long images through the 
B, R and I filters. Total exposure time amounts to 60 min 
in B, 75 min in R and 15 min in I. FWHM seeing was quite 
constant throughout the run with a mean value of 1.0". 
Images were corrected for bias and flat-fielded using stan- 
dard MIDAS procedures. All nights were of photometric 
quality. Observation of standard stars in M92 (Christian 
et al. 1985) allowed to calibrate the images with respect 



to the Kron-Cousins BVRI system. 

We also obtained low resolution spectroscopy of the 4 
brightest optical objects named A,B,C and D and located 
inside or close to the ROSAT all-sky survey error circle. 
Wavelength calibration was derived from the observation 
of Hg/Ar arc spectra and the instrumental response was 
computed using the flux standard star Feige 34. The Mul- 
tiple Object Spectroscopy mode of the OSIS instrument 
allowed the simultaneous acquisition of spectra from the 
4 objects through 1" slits and using the V150 grism. This 
configuration yields a FWHM resolution of 2nm and 
a useful spectral range of 365 to 990 nm. We acquired 6 
individual frames with exposure times ranging from 30 to 
60 min. The total spectral exposure time is 225 min. Spec- 
tra were flat-fielded, wavelength calibrated, extracted and 
corrected for atmospheric absorption and instrumental re- 
sponse using standard MIDAS procedures. 



Star 



B 



R 



I 



R- I 



A 
B 
C 
D 



23.83 ±0.20 



22.46+0.04 21.51+0.07 0.85+0.08 

22.53+0.04 21.03+0.05 1.32+0.06 

23.25+0.05 22.03+0.07 1.07+0.08 

21.90+0.03 21.34+0.05 0.52+0.06 



Additional images were obtained with the Keck LRIS 
(Low-Resolution Imaging Spectrograph) on 1999 February 
23. These observations took place under good sky condi- 
tions with a seeing around 0.9 arcsec. Two B and two R 
images of the field were taken with a total exposure time 
of 15 min in each filter. Using observations of one Lan- 
dolt standard star field, a photometric calibration with 
estimated uncertainty of 0.1-0.2 mag was achieved. 

4-. 2. Imaging data 

We show on Fig. Q the summed CFHT R band image with 
the 90% confidence level ROSAT HRI and ROSAT survey 
error circles over-plotted. The Keck B image is displayed 
in Fig. ||. We calibrated astrometrically our CCD image 
using 5 USN0-A2 stellar hke objects. The attitude cor- 
rected HRI 90% confidence error radius of 2" shown here 
includes the additional astrometric error of 0.7" arising 
from the CCD calibration. The only detectable object in 
the ROSAT HRI error circle is C. 

R band image profile measurements reveal that objects 
A and B are resolved. A seems extended in all directions 
whereas B has a stellar-like core with diffuse emission to- 
wards the SE direction. On the other hand, objects C and 
D appear unresolved. 

BRI photometry of objects A,B,C and D as derived 
from CFHT observations is listed in Table ^. Only object 
D is bright enough to be detected in the CFHT B band 
image. We estimate limiting magnitudes of B ^ 24.6 and 
R ~ 25.0 on the summed images. 

Keck photometry of object C gives R = 23.3 and B-R 
= 2.6 consistent with CFHT observations. The limiting 
magnitudes of the Keck images are estimated to be B ^ 
27 and R ~ 26. 

4-3. Spectroscopic data 

In general, the signal to noise ratio of the average spec- 
tra is not good enough to unambiguously measure a red- 
shift or detect the presence of weak emission lines. Telluric 
absorption lines (e.g. the A6800-6900A complex) are not 
pronounced. 

Identifying the flux drop bluewards of 6OOOA in object 
A (see Fig. ^) with the Ca break leads to a redshift of 
~ 0.5. At this redshift, other spectral features such as G 
band, H/3 and Mg band may be seen in the spectrum. 
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Fig. 4. The summed CFHT R band image showing the 
position of the ROSAT HRI attitude corrected error cir- 
cle (rgo = 2"), the ROSAT PSPC survey error circle (rgo 
= 15") and the 4 objects for which we obtained low reso- 
lution spectroscopy. 



The spectrum of object B leaves hardly any doubt 
that the stellar-like core is in fact a dwarf M3-4 star, the 
extended emission being then a likely background field 
galaxy. We show on Fig. |^ the observed spectrum together 
with that of a comparison M3V star extracted from the at- 



las of Torres-Dodgen & Weaver ( |1993| ). The Nal hnc, TiO 
and CaH bands are clearly detected. Neglecting the con- 
tribution of the background galaxy, the R-I colour index 
is also consistent with that of a M3-4V star. 

As object C is the only one detected in the small atti- 
tude corrected HRI error circle, special attention was given 
to its analysis. Being about 0.7 mag fainter than star B 
the spectral features of C are obviously less recognizable. 
There is however good evidence that C is also a M star, 
of slightly earlier spectral type than star B. For compari- 
son, we show on Fig. ^ the flux calibrated spectrum of C 
together with that of a template MOV star extracte d from 
the spectral atlas of Jacobi Hunter and Christian (1984). 
The NaD line and the broad TiO molecular bands visible 
in the smoothed MOV spectrum can be also seen in ob- 
ject C. The similarity of the energy distributions is also 
striking and the R-I and B-R colour indexes of C are also 
consistent with a ~ M2V star. All these evidences support 
the conclusion that C is a rather early M type star. 

Finally, based on Ca break, H/3 and NaD line positions, 
the spectrum of object D (Fig. ^ suggests an identification 
with a galaxy at z ~ 0.3. 



Fig. 5. The summed Keck B band image with the ROSAT 
PSPC survey error circle (rgo — 15") overplotted. North 
is at top and East at left. 
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Fig. 6. Flux calibrated spectra of objects A and D. Spec- 
tra have been smoothed by a Gaussian filter with a — 
ISA 



5. Discussion 

5.1. The nature of RX J1605. 3+3249 

The position of object C in the attitude corrected HRI 
error circle is somewhat puzzling but may not be fully 
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Fig. 7. Flux calibrated spectrum of object "B" (bottom) 
plotted together with that of a template M3V star (top) 
. The spectrum of B is smoothed using a Gaussian filter 
with cr = f SA 
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Fig. 8. Flux calibrated spectrum of object "C" (bottom) 
plotted together with that of a MOV (top). Both spectra 
were smoothed using similar Gaussian filters 



significant. At I = 53°, b = 48°, the surface density of 
stars brighter than R = 23.5 is 12,300 deg'^ (Robin & 



Creze 1986). There is therefore more than 1% probability 
that a star falls by chance in the 2" radius error circle. 

The photometric distance to the M2V star is ~ 7kpc. 
If the X-ray source were physically associated with the 
star, its radius would be i? ~ 80 km, i.e. too large for a 
neutron star and too small for a white dwarf. The emit- 
ting area could be compatible with a polar cap heated by 
accretion as seen in AM Her systems in which case, our op- 
tical observations could have been obtained during a low 
state. This explanation cannot be strictly ruled out, as the 
white dwarf could be cool enough in the low state {T^s < 
20,000 K) to remain undetected at B = 25.85 which is 
the B magnitude of object C. We also do not have enough 
spectral resolution to detect the emission lines which could 
reveal heating of the late M star by the white dwarf. How- 



ever, the absence of X-ray and optical variability and the 
unusually hot temperature of RX J1605. 3-1-3249 compared 
to those of black body components in polars make this 
possibility rather unlikely. 

A hot white dwarf would have to be located at unrealis- 
ticly large distances {d ^ 270 kpc), although this distance 
may be overestimated by the black body fit. The black 
body temperature is also hot compared to that observed 
from the hottest known PG 1159 stars (e.g. Werner et al. 
1996). Finally, if the high temperature of the white dwarf 
were due to nuclear burning at its surface, we would de- 
tect the heated accretion disc and mass donor star or the 
surrounding nebula such as in SMC N67. 

In general, all classes of soft emitters else than neu- 
tron stars are difficult to reconcile with the observational 
picture. Taking the expected V magnitude of the M2V 
star (V = 24.24) as an upper limit to the optical emission 
from the ROSAT source implies Fx/Fopt > lO''. Com- 
paring with the Fx/Fopt distribution of bright ROSAT 
all-sk y sur vey sources identified in SIMBAD (e.g. Motch 
et al. 1998) shows that all galactic or extragalactic classes 
of sources other than isolated neutron star s are probably 
ruled out (see also Fig. 3 in Schwope et al. 199!: ). In par- 
ticular, the source is optically too faint to be identified 
with even extreme cases of cataclysmic variables or AGN. 

This conclusion is independent of the HRI attitude cor- 
rection since none of the other optical candidates studied 
in the large ROSAT survey error circle is a likely counter- 
part of the X-ray source. In particular, the extragalactic 
objects lack the broad emission lines usually seen in soft 
AGN (e.g. Greiner et al. |1996D. 



Assuming a neutron star radius of 10 km implies a 
source distance of 900 pc for full surface emission. How- 
ever, this distance estimate is very sensitive to the actual 
energy distribution in the ROSAT band. For instance, Ra- 
jagopal and Romani (199^) have shown that blackbody 
fits to neutron star model atmospheres folded through the 
ROSAT PSPC tend to overestimate the effective temper- 
ature by a factor of up to 3 depending on chemical com- 
position. Applying the maximum correction could bring 
the source to much closer distances {d 100 pc). 

The total galactic Nh in 
RXJ1605.3+3249 is 2.47 lO^"' 
Lockman, 



1990) 



IS 

about 



twice 



the direction of 
' cm^^ (Dickey, J. & 
that derived from the 
blackbody fits to PSPC data (see Fig. However, 
this difference may not be significant since similar to 
the effective temperature, the estimated photelectric 
absorption sensitively depends upon the assumed soft 
X-ray energy distribution. 

We conclude that RX J1605.3-h3249 exhibits aU the 
features, soft thermal- like spectrum, high Fx/Fopt, ex- 
pected from an isolated neutron star. 



C. Motch et al.: The isolated neutron star candidate RX J1605. 3+3249 



7 



5.2. Accretion from the interstellar medium 

As for RX J0720. 4-3125, the average particle density avail- 
able for accretion in the vicinity of the neutron star is 
probably small. IRAS maps do not reveal any particular 
density enhancement in the direction of RX J1605. 3-1-3249. 
Using the mean particl e den sity versus scale height law 
of Dickey & Lockman (199C) yields average densities of 
only 0.3 cm^^ at 100 pc and 0.06 cm^'^ at 300 pc. In order 
to explain the X-ray luminosity of the source by Bondi 
Hoyle accretion, very low values of the total relative plus 
sound speed velocities must be achieved, 12 km s^"'^ and 
3.5 km s"""^ for d = 100 and 300 pc respectively. In partic- 
ular, large distances of the order of 900 pc as suggested 
by PSPC spectral fitting seem incompatible with the ac- 
cretion model. However, this major difficulty could be by- 
passed assuming the presence of a local overdensity or if 
the blackbody modelling strongly overestimates the total 
accretion luminosity. 

So far, three isolated neutron star candidates have 
been detect ed in the galactic plane (|&| < 20° (Haberl 
et al. |1998D and only two RX J1 605.3- t-3249 and IRXS 
J130848.6-h212708 (Schwope et al. |1999| ) above the plane. 
RXJ1605. 3-1-3249 is the X-ray brightest high galactic lati- 
tude candidate. The steep decrease of particle density with 
distance for RXJ1605. 3-1-3249 could allow to sensitively 
test the mechanism leading to X-ray emission in this newly 
discovered class of objects. 



5.3. Cooling neutron star 

The strongest argument in favour of a young cooling neu- 
tron star is probably the stability of X-ray emission over 
long time scales. This property of RX J1605. 3-1-3249 is also 
shared by RX J0720.4-3125 and RXJ1856. 5-3754. Neu- 
tron stars accreting in binaries always exhibit some kind 
of variability on a large range of time scales from seconds 
to months. Although the accretion conditions prevailing 
in binaries are usually far from those assumed for accre- 
tion from interstellar medium, it is still puzzling that none 
of the members of this new class so far studied in details 
shows convincing evidences for variability. 

For normal pulsars, the blackbody temperature of 1.1 
10^ K implies cooling ages in the range of 2x 10^ to 10^ yr 
depending on the presence of accreted material at the sur- 



face of the neutron star (Chabrier et al. 1997). Considering 
the probable overestimation of Tcff by the blackbody fit, 
ages of up to 10^ yr are still possible. 

RX J1605.3-h3249 and its cousins RX J0720.4-3125 and 
RX J1856. 5-3754 have in common two properties which 
are never encountered together in other classes of neutron 
stars, i) absence of strong radio emission and ii) absence 
of luminous hard X-ray tail above the thermal spectrum. 
These features could be used to define the new class of 
objects. 



Simultaneous absence of radio and hard X-ray emis- 
sion is understandable in the framework of old accreting 
neutron stars. The low magnetic field or the long spin pe- 
riod necessary for accretion to take place are likely to put 
the pulsar beyond the death line, in the graveyard, and the 
emitted X-ray spectrum is expected to resemble that of a 
soft black body for a large range of parameters (Zampieri 
et al. 1995). We envisage below the implications of such 



properties for young cooling INS. 
5.3.1. Radio emission 

The most simple explanation for the absence of radio emis- 
sion is that the radio beam does not cross the earth. The 
beaming fraction which is the proportion of the sky swept 
by the radio beam decreases with increasing spin period 
(Biggs 1990 ) and is of the order o f 0.2 for the overall pul- 
sar population (e.g. Lyne et al. 1998 ). Therefore, radio 
detection means may miss a large fraction of the pulsar 
population among which a part may be hot enough to be 
detected in soft X-rays. 

As noted by Kulkarni & van Kerkwijk (1998), 



this explanation cannot hold for the pulsating source 
RXJ0720. 4-3125 because the time needed to brake the 
neutron star to the long spin period of 8.39 s (assuming a 
magnetic field of lO^'^ G) is much larger than the cooling 
time. For the same reason, the absence of radio emission 
from RX J1605. 3-1-3249 is unlikely to be due to a position 
beyond the death line as this would also imply rather long 
spin periods incompatible with the hot emission. 

Another possibility is that RX J1605. 3-1-3249 is a mag- 
netar with a dipolar surface field B larger than about 
4xlO^'^G in whic h case radio emission may be quenched 
(Heyl & Kulkarni |l99|). 



5.3.2. X-ray spectrum 

The fact that the X-ray spectrum is to the accuracy of 
the measurements thermal-like suggests a much reduced 
magnetospheric activity compared to other known neu- 
tron stars. 

Among the 27 pulsars detected in X-rays and listed in 
Becker & Triimper (1997), only three middle aged pulsars 



(Geminga, PSR B0656-f 14 and PSR B1055-52) have in 
addition to a power law a recognizable black body com- 
ponent in their soft X-ray energy distribution. An interest- 
ing case is PSR B0656-t-14, a radio emitting pulsar about 
~ 10^ yr old with a magnetic field of 4.7x10^^ G and lo- 
cated at a distance of 760 pc. PSR B0656+14 exhibits a 
ROS AT PSPC spectrum (Tbb = 80-90 eV; Possenti et al. 



19961 ) strikingly similar to that of RX J1605.3-h3249. The 



additional faint hard component needed to fit the spec- 
trum of PSR B0656-I-14 would not have been detected in 
RXJ1605. 3-1-3249 because of the lower statistics. 

Neutron stars born with magnetic field B > 10^"* G, 
the magnetars, are thought to be powerful soft X-ray emit- 
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ters because magnetic field decay provides an additional 
source of heat (Thompson & Duncan 1996, Heyl & Kulka- 
The young magnetars associated with soft 7- 



1998) 



ray repeaters such as SGR 1806-20 or SGR 1900+14 ex- 
hibit powerlaw-like quiescent X-ray spectra without ev- 
idences for black body components. These non-thermal 
energy distributions could be the signature of a compact 



synchrotron nebula (Marsden et al. |1998). It has been pro 



posed that the class of anomalous braking X-ray pulsars 
coul d be related to magnetars (e.g. Thompson & Dun- 
can 1996) and could constitute a later, less active soft 



7-ray repeater phase. In general, anomalous X-ray pulsars 
have again powerlaw-like spectra with in two cases possi- 



ble blackbody components (see Thompson & Duncan |1996 



and references therein). Pulsating sources like RX J0720.4- 
3125 could represent an even later stage of magnetar evo- 
lution with remaining very high magnetic field (B ^ 10^^ 
G, Heyl & Hernquist 1 



and perhaps still the possi- 
bility to emit powerful 7-ray bursts on occasion. Because 
of the additional energy source a magnetar could reach 
the temperature of 1.1 10^ K after 10^ yr (Heyl & Hern- 
quist 1998| ). It is however unclear whether the absence of 
a strong non-thermal component in the X-ray spectrum is 
compatible with the remaining high magnetic field. 



5.4- Expected brightness of the optical counterpart 

Taking into account the unfortunate possibility of a chance 
alignement between the neutron star and object C implies 
a B magnitude fainter than - 26 for RX J1605.3+3249. 
Because of the relatively high temperature, the extrapo- 
lation of the black body seen in soft X-rays to the optical 
regime would imply extremely faint optical magnitudes 
close to V = 30. However, all neutron stars observed so 
far display optical continuum above the Rayleigh- Jeans 
tail of the soft X-ray thermal component. This is not unex- 
pected since black body fits tend to overestimate Toff • Fur- 
thermore, a non thermal optical component has been de- 
tected in at least two cases, Geminga and PSR B0656+14. 
Scaling the optical flux with the PSPC count rate of 
RX JG72G.4-3125 and RX J1856.5-3754 and neglecting any 
temperature effects yields V--27.2 for RX J1605.3+3249. 
On the other hand, if RX J16G5. 3+3249 is similar to 
Geminga or PSR B0656+14 its B magnitude could be as 
bright as our limit of 26. Therefore the source may well 
be bright enough to be optically identified with current 
means and optical imaging could allow the detection of 
proper motion which is a crucial test for determining the 
X-ray powering mechanism. 



source is due to thermal emission from a nearby isolated 
neutron star. However, based on the presently available 
data, it is not possible to distinguish between accretion 
from the interstellar medium or cooling as the main X-ray 
emitting mechanism. 

The constancy of the X-ray flux on various time scales 
and the difficulties encountered by the accretion model 
as a result of the small mean ambient densities are argu- 
ments, although not fully compelling, for a cooling neutron 
star. 

The undetermined neutron star spin period and the 
lack of sensitive measurement of a hard X-ray component 
prevent us from drawing any firm conclusion on the nature 
of the source. One possibility is that RXJ1605. 3+3249 is 
a twin of PSR BG656+14 but that the radio beam does 
not intercept the earth. Alternatively, RXJ1605. 3+3249 
could be a magnetar, maybe similar to RXJ0720. 4-3125. 

Further sensitive optical and X-ray observations with 
the XMM and AXAF satellites could help to unveil the 
real nature of this object. 
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